home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 2002 November / SGI IRIX Base Documentation 2002 November.iso / usr / share / catman / p_man / cat3 / SCSL / claed0.z / claed0
Encoding:
Text File  |  2002-10-03  |  4.5 KB  |  133 lines
  1.  
  2.  
  3.  
  4. CCCCLLLLAAAAEEEEDDDD0000((((3333SSSS))))                                                          CCCCLLLLAAAAEEEEDDDD0000((((3333SSSS))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      CLAED0 - the divide and conquer method, CLAED0 computes all eigenvalues
  10.      of a symmetric tridiagonal matrix which is one diagonal block of those
  11.      from reducing a dense or band Hermitian matrix and corresponding
  12.      eigenvectors of the dense or band matrix
  13.  
  14. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  15.      SUBROUTINE CLAED0( QSIZ, N, D, E, Q, LDQ, QSTORE, LDQS, RWORK, IWORK,
  16.                         INFO )
  17.  
  18.          INTEGER        INFO, LDQ, LDQS, N, QSIZ
  19.  
  20.          INTEGER        IWORK( * )
  21.  
  22.          REAL           D( * ), E( * ), RWORK( * )
  23.  
  24.          COMPLEX        Q( LDQ, * ), QSTORE( LDQS, * )
  25.  
  26. IIIIMMMMPPPPLLLLEEEEMMMMEEEENNNNTTTTAAAATTTTIIIIOOOONNNN
  27.      These routines are part of the SCSL Scientific Library and can be loaded
  28.      using either the -lscs or the -lscs_mp option.  The -lscs_mp option
  29.      directs the linker to use the multi-processor version of the library.
  30.  
  31.      When linking to SCSL with -lscs or -lscs_mp, the default integer size is
  32.      4 bytes (32 bits). Another version of SCSL is available in which integers
  33.      are 8 bytes (64 bits).  This version allows the user access to larger
  34.      memory sizes and helps when porting legacy Cray codes.  It can be loaded
  35.      by using the -lscs_i8 option or the -lscs_i8_mp option. A program may use
  36.      only one of the two versions; 4-byte integer and 8-byte integer library
  37.      calls cannot be mixed.
  38.  
  39. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  40.      Using the divide and conquer method, CLAED0 computes all eigenvalues of a
  41.      symmetric tridiagonal matrix which is one diagonal block of those from
  42.      reducing a dense or band Hermitian matrix and corresponding eigenvectors
  43.      of the dense or band matrix.
  44.  
  45. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  46.      QSIZ   (input) INTEGER
  47.             The dimension of the unitary matrix used to reduce the full matrix
  48.             to tridiagonal form.  QSIZ >= N if ICOMPQ = 1.
  49.  
  50.      N      (input) INTEGER
  51.             The dimension of the symmetric tridiagonal matrix.  N >= 0.
  52.  
  53.      D      (input/output) REAL array, dimension (N)
  54.             On entry, the diagonal elements of the tridiagonal matrix.  On
  55.             exit, the eigenvalues in ascending order.
  56.  
  57.      E      (input/output) REAL array, dimension (N-1)
  58.             On entry, the off-diagonal elements of the tridiagonal matrix.  On
  59.             exit, E has been destroyed.
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. CCCCLLLLAAAAEEEEDDDD0000((((3333SSSS))))                                                          CCCCLLLLAAAAEEEEDDDD0000((((3333SSSS))))
  71.  
  72.  
  73.  
  74.      Q      (input/output) COMPLEX array, dimension (LDQ,N)
  75.             On entry, Q must contain an QSIZ x N matrix whose columns
  76.             unitarily orthonormal. It is a part of the unitary matrix that
  77.             reduces the full dense Hermitian matrix to a (reducible) symmetric
  78.             tridiagonal matrix.
  79.  
  80.      LDQ    (input) INTEGER
  81.             The leading dimension of the array Q.  LDQ >= max(1,N).
  82.  
  83.      IWORK  (workspace) INTEGER array,
  84.             the dimension of IWORK must be at least 6 + 6*N + 5*N*lg N ( lg( N
  85.             ) = smallest integer k such that 2^k >= N )
  86.  
  87.      RWORK  (workspace) REAL array,
  88.             dimension (1 + 3*N + 2*N*lg N + 3*N**2) ( lg( N ) = smallest
  89.             integer k such that 2^k >= N )
  90.  
  91.             QSTORE (workspace) COMPLEX array, dimension (LDQS, N) Used to
  92.             store parts of the eigenvector matrix when the updating matrix
  93.             multiplies take place.
  94.  
  95.      LDQS   (input) INTEGER
  96.             The leading dimension of the array QSTORE.  LDQS >= max(1,N).
  97.  
  98.      INFO   (output) INTEGER
  99.             = 0:  successful exit.
  100.             < 0:  if INFO = -i, the i-th argument had an illegal value.
  101.             > 0:  The algorithm failed to compute an eigenvalue while working
  102.             on the submatrix lying in rows and columns INFO/(N+1) through
  103.             mod(INFO,N+1).
  104.  
  105. SSSSEEEEEEEE AAAALLLLSSSSOOOO
  106.      INTRO_LAPACK(3S), INTRO_SCSL(3S)
  107.  
  108.      This man page is available only online.
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.